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1.0 INTRODUCTION 

There Is a continuing interest in the determination of the density funotion of 
a material body. The usual approach is to attempt to r, ? oover the density funo- 
tion from knowledge of the total mass and the Newtonian potential (refs. 1, 2, 
3, and 4). It is well known that, in general, one cannot determine the den- 
sity function from knowledge of only these two items. 

In „hia work we postulate the existence of higher-order potentials, and demon- 
strate that the density function can be determined from this set of potentials. 
Moreover, we shall show that there is a one-to-one correspondence between such 
potentials and the density funotion, provided it is continuous. 

Finally, we note that this demonstration includes a method of constructing the 
density function from the set of potentials as well as indicating how one might 
approximate the density from incomplete knowledge of the potentials. 


2.0 MODEL AND NOTATION 


Suppose A is a material body. We model this by assuming that there is a 
nonnegative continuous function P defined on Euclidean three space denoted by 
E such that P has compact support K, where K has the same geometric shape 
as A and P is the density function for A. This notion precludes A 
having "point masses," but does not require that A be "in one piece;" i.e., 

K need not be connected. 

We denote the real line by R, and r,he norm on E, induced by the inner product 
<•,*>, by |M|. 

If r is a positive number and P is a point in E, then 
S P (P) = iQ in E: | |p - q| | < r) 


and 


3s r (P) = (Q in E: | |p - q| | = r) 


3-0 MATHEMATICAL ANALYSIS 

Suppose K is a compact set in E, which is the support for a nonnegative con- 
tinuous function P, which is defined on E. Let C be a point in E and r 
be a positive number such that K is a subset of S r (C). Let D denote the com- 
plement of S r (C). 
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Definition 1. The statement that M is the mass of K means M is the 
number 


J p(Q)dQ 

E 


We understand that the integral is a triple integral. 

Definition 2. For eaoh nonnegative integer r, V n is the real valued function 
defined on D by 


f P(Q) 

V n (P) = J dQ for P in D 

S r (C) 1 1 P - Q j | n 

Notice that for each P in D, Vq(P) = M, and also that V-|(P) is the Newtonian 
potential at P due to p, 

For each positive integer n, we shall refer to V n as the nth-order potential 
and the set { V n : n = 2, 3» •••} as the higher-order potentials. 

We now state the main purpose of this work. 

Theorem. Suppose p is a nonnegative continuous function defined on E, which 
has compact support K. Si is a ball containing K, and Sg is a ball 
properly containing Si and is concentric with There is a one-to-one comes- 

CO 

pondence between the set V = { {V n (P)} n= Q. P is in the boundary of a hemi- 
sphere of Sg} and p; moreover, p can be constructed from V. 

Definition 3. Suppose p is in D. Let mp be the function defined on R 
by 


rop(x) = 0 


if x < 0 

mp(x) = J* 

p(Q)dQ 

i^ x > 0 


S X (P) 


2 


\ 
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nip(x) s / p(Q)dQ if x > 0 
S X (P) 


Figure 1.- The function mp. 


Observe that for each point P in D, mp is a real-valued nondecreasing 
continuous function on R, which has a continuous derivative given by 


1 

mp(x) 

= 0 

if x < 0 

mp(x) 

= J p(x,z)dz 

if S r (C) intersects 3S X (P) (S r (C)nS x (P) ) 

s r (c)nas x (p) 


f 

mp(x) 

= 0 

if x > 0 and 3S X (P) does not inter- 
sect S r (C) 

Definition 4. 

If P is a point in D, 

then 


up = inf .{ | J P - Q| j : Q is in S r (C)} 

vp = sup . { | | P - Q( | J Q is in S r (C)} 

Notice that up and vp are readily determined for each P in D. 
Using definitions 3 and 4, observe that 


3 
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VP 

v n (P) * f — dmp(t) for n = 0, 1, 2, ... 
t« 

Up 

Making a change of variables s = 1/t for up < fc < vp> we have 
vp-1 

V n (P) = / s^dinpd/s) for n = 0, 1, 2, ... 

up-1 

Definition 5. For each point P in D 

A P = vp -1 
Bp = up -1 


and 


Mp(x) = M - mp(1/x) 


for Ap < x < Bp 


Observe that Mp(Ap) = 0, and also that Mp(Bp) = M; hence, we may extend Mp 
to (0, Bp) by defining Mp(x) = 0 for 0 < x < Ap. Note also that 

. » 2 , 
Mp(Ap) = Mp( 1/xp) ( 1/Ap) = 0; hence, Mp has a continuous derivative on (0, 

B P ). 


We now have 


B P 

V n (P) = J x n dM P (x) for n = 0, 1, 2, ... 
0 

» i 2 

Notice also that Mp(Bp) = Mp( 1/Bp) ( 1/Bp) = 0. 
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We seek to arrange matters such that we are integrating over the interval 
(0,1). We have two options at this point. If we restrict our attention to 
only those points P outside of Sp+i(C), then Bp * up" 1 <1. We then are 
able to extend Mp to (0, l) by setting Mp(x) a M, if Bp < x < 1, ard to 
observe that Mp has a continuous derivative on (0, l), and also~that 


1 

V n CP) * / x n dMp(x) for n * 0, 1,2, *** 
0 


If we do not elect to make this restriction, then by letting y * x/Bp for 
0 < x < Bp, we have 


1 

V n (P) s Bp / y n dMp(Bpy) for n = 0, 1, 2, **• 
0 


If we now define V n (P) s Bp" n V n (P) for P in D and n = 0, 1, 2, and 

define Mp(x) = Mp(Bpx) for 0 < x < 1, then we have 


1 

V n (P) = J x n dMp(x) for n = 0, 1, 2, ••• 
0 


let us proceed under the assumption that Bp < 1 , keeping in mind that Mp so 
defined is a nonnegative nondeereasing function that has a continuous deriva- 
tive. 


00 

If P is in D, we can construct a sequence of functions {M n } n --| from 

co oo 

the number sequence {V n (P)}n=;0 suc *' 1 that fche function sequence {M n } n= i con- 
verges uniformly to Mp on (0, l), and there is a one-to-one correspondence 

OQ 

between Mp and { V n (P ) > n -o * (This construction is further defined in appen- 

( 00 

dix A.) Hence, we have a one-to-one correspondence between Mp and { V n ( P ) } n _ 
for each P in D. 


Let h > r + 1 and let us restrict our attention to only those points P 
that are in B, the boundary of a hemisphere of S^CC). For each point P 

let Lp denote the line containing P and C, 

Definition 6. For each number x where S X (P) interseots S r (C) f let 

wp(Q), for each point Q in gS x (P)r)S (C), be the point on Lp common to 

S r (C) and 3S X (P). 




Me) w'q) 


Figure 2.- The function wp(Q), 


Definition 7. For each point P in B, and each number x, let 


G(P , x) s mp(x) 


Lemma. G(P,x) is continuous and bounded on Bx (-oo, co). 

The proof is given in appendix B. Notice that for P in B, ||P 
Definition 8. For each point P in B, let 


Hp(x) = mp(x + h) 


for x in ( -oo, oo) 


With this definition we observe that Hp(0) = mp(h), and also that 
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00 

/ Hp(x)e“i x ^dx 


-00 


s / p(Q)e-i<wp(Q).y>dQ for y in L P 
E 


where d is the number suoh that y = du with u being a fixed unit vector 
at C, along Lp. 


The first of the above integrals is (2ti) 1/2 Hp(d) , where H P is the Fourier 
transform of Hp. 

Since H P is zero (except on, at most, a hounded set) and Hp is continuous, 

A 

it follows that H p is a continuous Li(-», w) function for each P in B 

y\ 

(refs. 5 and 6). Moreover, if follows that H P (d), when viewed as a two- 
place function, is bounded and continuous and is in L-](Bx(-», <®)), as shown 
by the lemma. 

For each point Y in L P , let 

pp(y) = (2 tt)-3/2 f p(Q) e ~i < wp(Q), Y>dQ 
E 


Definition 8. For each point P in B, let T P be the transformation from E 
onto E, which is defined as follows. For each number x let tt x be the plane 
orthogonal to Lp, containing the point D x = (1 - x/||P - Q||)P +• (x/ 1 j P - Q| | )C. 
If D x is not in S r (C), then Tp is the identity map on D x . If D x is in 
S r (C), then for each point Q in S r (C)r)3S x (P) , T P (Q) is the orthogonal pro- 
jection of Q onto tt x . If Q is in the plane containing 3S r (C)r|8S x (P) and 
Q is not In S r (C), then TpQ is the orthogonal projection of Q onto tt x . 


7 


Figure 3s" The transformation Tp. 


Notice that for each P in B, Tp is a reversibly continuous mapping from E 
onto E, which leaves each point in Lp fixed. 

If Q is in S r (C) n aS x (P) and Y is in Lp, then 


<.TQ,TY> = <TQ,Y> 

» ll TQ ll 11*11 9 

= ||W|| oo» 0 ||Y|| 

= ] | wp<Q> | | ||Y|| 

= <wp(Q) , Y> 


We can now write 



p P (Y) » (2ti)“3/2 j p (TQ)e“i< T, 5|Y>dTQ 
TE*E 

s (2tt)"3/2 j p(Z)e“ i<z » Y >dZ 
E 

where Z = TQ and V is in Lp. 

Hence, 

pp(Y) s pp(Y) 

the Fourier transform of p for Y in Lp. 

Therefore, 

pp(Y) = ( 2 tt ) “ 1 H P (d) for Y in Lp 
where d is the number such that dU = Y. 

Notice that if Y i C, there is only one point P in B such that Y is 
in Lp. 

A 

Definition 9. Let H be the function defined on E by 
H(Y) = Hp( Y) 

if Y A C and Y is in Lp and H(C) =0. 

A 

Definition 10. Let p be the function defined on E by 
p(Y) = (27r)-lH(Y) 
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In review, if for each point 

f 

mine uniquely mp(x) for x 
H(Y) for all Y in E and, 


00 

P In B we know {V n (P)} ns o> then we can deter 
in (-», °° ) . From this we are able to determine 

A 

hence, p for all Y in E. 


Continuing, suppose Q is in S r (C) and Q is not C. Then there is a 
unique point P and B such that Q is in the line Lp and 


p(Q) e (2 t t)~3/2 f p(R)e i<w P^)»Q>dR 
E 


Recall that p = (2 tt)“^ H from which we have 


p(Q) 


(2tt)“5/2 


/ H(R)e i<w P(R)>Q>dR 
E 


Hence, for each ’ Q in S r (C), Q 4 C, we can determine p(Q) uniquely, and 
therefore, since p is assumed continuous, we have determined p uniquely and 

00 

the argument is complete. Hence, if {V n (P)} 11= o is known for each point P 
in B / the p is uniquely determined under the assumptions that p is 
continuous. 


1) . C COMMENTS 


The preceding analysis could have been carried out similarly had we elected 
not to assume that Bp < 1 , 

Also, the requirement that P belongs to B, the boundary of a ball containing 
K, is not necessary, What is essential is that for each point Q in S r (C) 

Q not C, there be a point P not in S r (C) such that Q is in the line Lp. 


5.0 APPROXIMATION OF DENSITY 

, N 

Suppose that for each point P in B, only {V n (P)} n _o is known. Then one 
can only construct, for each P in B, the first N approximates to Mp. 

This then leads to only the first N approximates to Mp and, hence, to Hp. 
One can still proceed, using only these first N approximates since each is 
bounded by a summable function. Notice that for P in D 

t 


10 


80FMl>2 


1 

V n (P) = / x n dMp(x) for n = 0, 1, 2, ••• 

0 

and, hence 

1 

V n (P) = ( 1 / ( n+ 1 ) ) / Mp(x)dx n+ 1 
0 

1 

= ( 1/ (n+1 ) )Mp(1 ) - (1/(n+1)) / x n +1dMp(x) 

0 

for n = 0, 1 , 2, • • • 

Therefore, if we let 

D n (P) = Mp( 1 ) - (n+1 ) V n (P) for n = 0, 1, 2, ••• 
and P in D, we could, using the construction in appendix A, generate a 

— 00 'i 

sequence {M n } n=0 of functions on (o, l), which will converge to Mp on 
(0, 1}. 

A second notion is that perhaps for only a finite number of points P in B 
00 

is {Vn(P))n=0 known. Then one could proceed as in the argument except that 

A 

one interpolates values for Hp for those points P in B for which <» 
{V n (P)} n=0 are not known; then proceed as before. 

Both of these notions, as well as the convergence, rates should be addressed. 
6.0 ADDITIONAL PHYSICAL NOTIONS 

Suppose that for each nonnegative integer n we define 
F n (P) = grad V n (P) for P in D 


11 
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Then Fq(P) s 0 and F i( p ) is fch 0 usual inverse square law. Notice also that 


F n (P) = —n / 
S r (C) 


><Q) 


PrQ 


|P-Q|| n+1 | |P-Q| | 


dQ 


for P in D and ns 1 , 2, • • • 

If g is a real valued analytic function on R, 


g(x) s i a n x n 
nsi 


then 


V g (P) - l a^P) 

nsi 


is well defined, as is 


F g (P) » Z a n F n (P) 
nsi 


This is readily observed by noting that 


N N 

S a„V„(P) = / p(0) E a n /||P-Q||» d8 

n= 1 S P (C) nsi 


and, hence 


V p > = , Ls p(Q) g<1/ l l p - Q l I )d ° 


for P in D. 


12 
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In a simlliar fashion 


N 

E a n p nt p ) 
n=1 



P-Q 

i l p " Q l I 


dQ 


f P (Q) 
S r (C) | |P-Q 1 1 2 


N 

E a n (-n) ( 1/ j jP-Q j j ) n_1 
n=1 


P-Q 

ll p -«l 


dQ 


and, henoe 


F g (P) 


f P (Q) 
sjc) | jP-Q | J 2 


g* (1/ | |P-Q | )) 


P-Q 

!l p - Q H 


dQ 


Let us denote Vg as the sum potential and Fg as the sum force, each rela- 
tive to g. Notice that if g(x) = ax, then we have the familiar Newtonian po- 
tential and force (ref. 7). 

An interesting question that should be addressed follows; if 


V(P) 


00 

E 

n=1 


a n V n (P) 


is well defined for all P in D, then is there an entire function g such 
that 


g(x) = E a n xn 
n=1 


13 
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7.0 CONCLUDING REMARKS 

We have demonstrated that if A is a material body (not necessarily having 
only one component), which has a continuous density function, and A is 
contained in a sphere S, then if one knows all the potentials (Newtonian and 
higher orders) for A at each point P of a sufficiently dense set B, then 
the density function for A can be determined uniquely. 

In addition, a notion of sum potential and sum force has been introduced. 

There is clearly a question as to how to determine the higher-order forces and 
potentials from experimental data as well as to determine the coefficients of 
g, if such a g exists (ref. 8). 

The notions mentioned here should not be considered limited to density func- 
tions of material bodies. If one chooses the alternative approach (i.e., not 
requiring that Bp < 1), then one could consider questions dealing with forces 
or fields near to sources of such phenomena. It should be noted that the domi- 
nant terms then appear to be the higher-order potentials and forces, whereas in 
dealing with questions about effects far from the souroe the higher-order poten 
tials and forces appear to play a far smaller role. The material in reference 
8 is very interesting with regard to the above comments. It would also be 
interesting to consider the work in reference 9 in regard to the notions 
introduced here, 
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APPENDIX A 

GENERATING A SEQUENCE OF FUNCTIONS 


OC 

We Indicate here a method of generating a sequence of functions {M n } n -i, 
each of which is real valued and defined on (0> 1), from a number sequence 

{ c n}n=0* (Please refer to references A-1, A-2, A-3, A-4, A-5, and A-6.) 

00 oo 

Given a number sequence {C n }n=0» the function sequence {Mn} n=1 is called 
the associated function sequence. 

Suppose n is a positive integer and x is a number in the number interval 
(0, 1]. There is a unique integer k(n,x) such that 


k(n,x)/n £ x < (k(n,x)+1)/n 


For each positive integer n define a step function M n on fO, 1) as 
follows: 


M n (0) = 0 


k(n,x) n n-t n _ t 

M n ( x ) = Z (t) 2 ( i )(-1) x C i+t 

t=0 L=0 


for x is in (0, 1). 

Reference A-4 shows that if M is a real valued function on (0, 1), which has 
a continuous derivative, and for each nonnegative integer n 


1 

C n = J x n dM(x) 
0 


00 

then the associated function sequence {M n } n= i converges uniformly to 
M -M(0) on (0, 1). 


A-1 
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APPENDIX B 

ESTABLISHING THE LEMMA 


The lemma la established here. For all P in B and x In (-»* 00 ) 
G(P,x) s Mp(x) 


If x < 0 G(P,x) = 0 for all P in B. 

If x > 0 and aS x (P) does not intersect S r (C), then G(P,x) = 0. 
Suppose now that S X (P) intersects SpCC) 


j G( P , x) 


J p(x, z) dz| 

S r (C)n3S x (P) 


and, hence 


< | p | / dz l 

S r (C) n 9S x (P) 


j G(P,x) | < |p| A 

for all x in (-°°, «>) and P in B. 

Let us now demonstrate continuity. If x is a number such that S r (C)n3S x (P) 
Is a singleton, then G(P,x) = 0. 

If 3S X (P) does not intersect S r (C), then 0(P,x) = 0; and if x < 0, then 
G(P,x) = 0. 

Suppose each of ? and P is in B, and each of x and x is a positive 
number. 


B-1 
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|G(P,x) - G(P,x)| s f p(x, z)dz 

S (C) n9 S (P) 
r x 


/ „ p(x, 

S (C) n gS (?) 
r x 


z)dz 


<, f |p(x, z) - p (x, z) I dz 

s (C) n3 s (P) 

r x 




i 


(C) n9 s (P) 

X 


p(x, 


z)dz 


- J p(x, z)dz 

s (c) n asjp) 

r x 


* 


p(x, 

(OnasjP) 
r x 


z)dz - f p(x, z)dz 

S (C) n gS__(P) 
r x 


Since p is continuous, we may make the first integral as small as we please 
by choosing |x - x| small. Since p is continuous and bounded, we may make 
the second terra as small as we please by choosing |x - xj sufficiently small 
with P and x fixed, and by satisfying each of the above we can make the 
last term as small as we please by choosing ||P - P|j sufficiently small. 


B-2 
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